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Power spectrum of HI intensity fluctuations in DDO 210 
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ABSTRACT 

We measure the power spectrum of HI intensity fluctuations in the extremely faint (Mb ~ — 
10.9) dwarf galaxy DDO 210 using a visibility based estimator that is well suited in very 
low signal to noise ratio regimes. DDO 210's HI power spectrum is well fit by a power law 
Pm(U) = AU a with a = -2.75 ± 0.45 over the length-scales 80 pc to 500 pc. We also 
find that the power spectrum does not change with an increase in the velocity channel width, 
indicating that the measured fluctuations correspond mainly to density fluctuations. However, 
Kolmogorov turbulence (with a velocity structure function spectral slope of 2/3) cannot be 
ruled out from the present observations. The value of the slope a is similar to that obtained in 
the Milky way. In contrast to the Milky way, DDO 210 has three orders of magnitude less HI, no 
spiral arms, and also no measurable ongoing star formation. The fact that the power spectrum 
slope is nonetheless similar in these two galaxies (and also similar to the values measured 
for the LMC and SMC) suggests that there is some universal, star formation independent, 
phenomenon responsible for producing fine scale structure in the gas. 
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^ 1 INTRODUCTION 



Evidence has been mounting in recent years that turbulence plays 
an important role in determining the physical conditions of the neu- 
tral ISM as well as for generating the hierarchy of structures seen in 
it ( see Elmegreen & Scalo (2004) and Scalo & Elmegreen (2004) 
for recent reviews). Observational evidence includes the fact that 
the fluctuation power spectrum of a variety of tracers (HI 21cm 
emission intensity, HI 21cm optical depth, dust emission) is a scale 
free power law. The slope of the power law of HI 21cm emission 
and absorption in our own galaxy, HI 21cm emission from the LMC 
and SMC are all ~ -3 (Crovisier & Dickey (1983); Green (1993); 
Deshpande et al. (2000); Elmegreen et al. (2001); Stanimirovic et 
al. (1999)). Similarly, the HI distribution in several dwarf galax- 
ies in the M81 group appears to be fractal (Westpfahl et al.(1999)), 
consistent with there being no preferred length scale. This scale free 
behaviour is characteristic of a turbulent medium; it is believed that 
the turbulence itself is driven by a combination energy input from 
spiral arms/bars and energy input from stellar sources (supernovae, 
stellar winds etc.). 

Previous quantitative measurements of the HI emission power 
spectrum have been limited to the few cases (the Milkyway, LMC 
and SMC) where the HI signal is extremely strong. In this paper, we 
discuss the power spectrum of the fluctuations of HI 2 1 cm emission 
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intensity in DDO 210, the faintest (Mb ~ —10.9) gas rich member 
of the local group. We use a visibility based estimator of the power 
spectrum; this is well suited to the current problem where the sig- 
nal is buried deep in the noise. Visibility based methods also have 
the advantage of having well understood statistics, and being free 
of the uncertainties involved in gridding and image deconvolution. 
The rest of the paper is divided as follows. The power spectrum es- 
timator that we use is discussed in detail in Sect. [2] the DDO 210 
data presented in Sect.|3| while the results of applying the estimator 
to this data are shown in Sect.|4]and discussed in Sect.|5| 



2 A VISIBILITY BASED POWER SPECTRUM 
ESTIMATOR 

The power spectrum of HI emission Phi (u, v) is the Fourier trans- 
form of the autocorrelation function £(l — I' , m — m') of the fluc- 
tuations in the HI brightness distribution 51(1, m) i.e. 

£(l-l',m- rri) = (51(1, m)5I(l' , m')) (1) 

and 

Pm(u,v)=J j £(l,m)e- l27T<ul+vm) dldm (2) 

Here (I, m) refers to directions on the sky and (u, v) to inverse an- 
gular separations. The Fourier relation in Eqn.[2]assumes that the 
angular extent of the galaxy is small, a spherical harmonic decom- 
position would be needed otherwise. Further, it is assumed that the 
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statistical properties of the small scale fluctuations in the HI distri- 
bution are homogeneous and isotropic. The angular brackets denote 
an average across different positions and directions in the galaxy. 
It follows that the HI power spectrum Pm(U) is a function of only 
the magnitude U = y/u 2 + v 2 . 

Now, the complex visibility function measured by an interfer- 
ometer is the Fourier transform of the brightness distribution, 



V{u,v) 



Til \ —i2-ir(ul4-vm) 17 1 

1(1, m) e K ' dl dm . 



(3) 



It should be noted that (u, v) in Eqn.[3]refers to the projected base- 
lines or antenna separations in units of the wavelength of observa- 
tion. The Fourier relation in Eqn. fallows us to associate an inverse 
angular scale with every baseline (u, v) which we also denote using 
a two dimensional vector U. Since the complex visibility measured 
by an radio interferometer is the Fourier transform of the HI inten- 
sity distribution, the squared modulus of the visibility is a direct 
estimator of the power spectrum, 



{V(U)V*(U))=Pm(U) 



(4) 



where (...) now denotes an average over all possible orientations of 
the baselines U. 

This estimator has been used by both Crovisier & 
Dickey(1983) and Green(1993) to measure the power spectrum of 
HI fluctuations in our own galaxy. Eqn.|5|is an idealised one, in 
practice the interferometer response differs from the Fourier trans- 
form of the sky brightness both because it is modified by a window 
function and because there is added noise. We first consider the ef- 
fect of the window function. In situations like in our own galaxy, 
where the HI emission fills the whole telescope beam, the window 
function is simply the primary beam of the telescope. In situations 
where the HI emission is from, for e.g., a distant galaxy whose an- 
gular extent is much smaller than the primary beam, the window 
function corresponds to the total extent of the HI emission. We will 
assume that, in general, the HI brightness distribution at any fre- 
quency v can be written as 



I v {l,m) = W v {l,m) [J„ + 6I„(l,m)], 



(5) 



where WvQjn) is the window function which characterises the 
large scale HI distribution in the galaxy at frequency v, /„ is the 
mean HI brightness and 5I v (l,m) characterises the small scale 
fluctuations. For a window function with angular extent 1/D, we 
expect its Fourier transform, W(U), to be sharply peaked at U = 
and to fall off rapidly for | U D. Then, at long baselines, where 
we can ignore the contribution from W v (I, m) I v , the square of the 
complex visibility is the power spectrum Phi(U) convolved with 
the square of the Fourier transform of the window function, i.e. 



{v v {u) v;{u)) = w;{u) ® Pbx{u) ■ 



(6) 



Both Crovisier & Dickey (1983) and Green(1993) ignored the win- 
dow function in their analysis, which is justified at large baselines 
if W„(U) is sharply peaked at U = and decays much faster than 
Pm(tf). 

We now look at the effect of measurement noise. The mea- 
sured visibility can be written as V V (U) = S V (U) + iV„, where 
S V (U) is the HI signal and N v is the noise. We then have 



(Vu(u) v;(u)) = {s v (u)s:(u)) + (n„n;) 



(7) 



where we see that squaring the complex visibility makes the noise 
bias positive-definite. For HI emission from our own galaxy this 
noise bias is small compared to the signal. On the other hand, 
the noise bias can be orders of magnitude larger than the desired 



power spectrum in the case of faint external galaxies. In principle, 
one could estimate (iV„iV^) from a line free frequency channel 
and subtract it from (V V (U) VJ(U)) for a channel with emission. 
However this requires that the noise statistics at different frequency 
channels be determined at an extremely high level of precision. Un- 
certainties in the bandpass response usually render such high levels 
of precision unachievable. 

The problem of noise bias can be avoided by correlating the 
visibilities at two different baselines for which the noise is expected 
to be uncorrected. The correlation between the visibilities mea- 
sured on two slightly different baselines is: 



Phi (U, AU) = {V V {U) V V *{U + AU)) 



Phi (U ) [W„(U - U ) X Wt(U + AU - U )] d 2 U , (8) 



Since W(U) falls off rapidly for | U |> D, it then follows 
that the two shifted window functions in Eqn. [8] will have a sub- 
stantial overlap only if \AU\ < D. Thus visibilities at two differ- 
ent baselines will be correlated only as long as \AU\ < D, and 
not beyond. In our analysis we restrict the difference in baselines 
to \AU\ « D so that W U (U + AU - U ) w Wu{U - U'). 
The expected value of the visibility correlation now no longer de- 
pends on AU allowing us to introduce the notation Pm{U) = 
(V V (U) V*0 + AU)smd 



Pm(U) = / P m (U ) | W V (U- 



U ) | 2 d 2 U , 



(9) 



We use the real part of the measured visibility correlation 
Phi (£0 as an estimator of the HI power spectrum Pm(U). This 
is the sum of the HI signal and a system noise contribution. The 
latter is unbiased and has a zero mean, its variance being deter- 
mined by the noise statistics of the individual visibilities. The imag- 
inary part of Phi (£0 is noise dominated. There will be a very 
small contribution from the HI signal because the assumption that 
W V (U + AU — ff) w W v {p -U)is not strictly valid. The re- 
quirement that the imaginary part of Phi (U) should be small pro- 
vides a self-consistency check and allows us to determine the range 
of validity of our formalism. 

A further simplification is possible for our estimator if we as- 
sume that W V (U) decays much faster than variations in Pm(U ). 
Such an assumption is justified at large baselines U 3> D if, for 
example, Phi (U) is a power law. We then have 



Pm(U) =Pm{U) 



W V (U 



(10) 



where the quantity in square brackets [..] in Eqn.^j|is a constant. 
We note that while Phi (U) directly estimates Phi at large baselines 
U 3> D, it is necessary to account for the convolution (Eqn.|9j at 
small baselines U ~ D. 

The utility of using visibility correlations as a statistical esti- 
mator of the HI power spectrum has been discussed previously by 
Bharadwaj & Sethi (2001) in the context of the large scale HI dis- 
tribution at high redshifts. Visibility based statistical techniques are 
also of interest in the efforts to detect the e poch of r eionization 
through interferometric HI obs ervations (eg. iMorales and Hewitt! 
l2004lBharadwai and Alil l2005 ) , and are also used in the analysis 
of inte rferometric observations of the CMBR (e.g. Hobs on et alJ 

mn). 
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3 DATA 

We apply the formalism developed above to interferometric obser- 
vations of DDO 210. Since DDO 210 is relatively close (at a dis- 
tance of 950±50 kpc Jlxe et all f 1999)) the observations have good 
spatial resolution as well as comparatively good signal-to-noise ra- 
tio (SNR). Further, DDO 210 has been observed with both the Gi- 
ant Metrewave Radio Telescope (GMRT) (Swarup et al. 1991) and 
in multiple configurations of the Very Large Array (VLA). A com- 
parison of the power spectrum derived from these two independent 
data sets allows a robust check against systematic biases. 

The GMRT data for DDO 210 have been discussed in detail 
in lBegum & ChengaluJ 120041) . We state the main observational re- 
sults b riefly below, the reader is referred to iBegum & ChengaluJ 
120041) for a more details. The HI emission from DDO 210 spanned 
the central 25 channels of the 128 channel spectral cube (with chan- 
nel width ~1. 7 kms -1 ). The HI disk of the galaxy is nearly face on. 
On large scales, the HI distribution is not axisymmetric; the inte- 
grated HI column-density contours are elongated towards the east 
and south. The maximum angular extent of the HI distribution for 
the galaxy is D ~ 0.8 kA and does not change significantly in the 
central 10 frequency channels. 

The VLA observations have been discussed in lYoung et alJ 

(2003). For the current analysis we use raw data downloaded from 
the VLA archive and reanalysed by us. 



4 RESULTS 

We estimated Phi(J/) by correlating all the baselines within 
| A[7 | ^ 0.2 kA. This value of A U was chosen as a compromise 
between having sufficient visibility pairs to get a good SNR and en- 
suring that AU itself remains small compared to D ~ 0.8 kA. All 
self correlations (i.e. correlations of a visibility with itself) were 
excluded. The measured correlations were then binned in bins of 
0.5 kA width. Each frequency channel was treated separately, and 
the analysis was also carried out on the line free frequency channels 
as a check. 

Fig-0A] shows the absolute value of the real and imaginary 
parts of Pm(U) for a channel with HI emission, while Fig. 0T4] 
shows the same quantities for a channel with no emission. For the 
channel with emission, we find that at baselines less than ~ 5 kA 



the real and imaginary parts of Phi (U) are both positive, and the 
real part is substantially larger than the imaginary part. For the 
channel with no emission the real and imaginary parts are small 
and of comparable magnitude. This is in agreement with the cross- 
checks discussed in Sec.|2|for determining the reality of the signal. 
Henceforth, we will ignore the imaginary part of Phi({7) and use 
the symbol Phi (U) to refer to only the real part. Note that the ap- 
parent increase in the value of Phi (U) seen at the longer baselines 
in Fig.0A]&[B] is because the noise is large in this range and 
we are plotting the absolute value of the estimator. The increase in 
noise is a consequence of the sparse sampling of the outer parts of 
the uv-plane. 

Fig.[2]shows Phi (U) determined separately from the GMRT 
and VLA data (for channels corresponding to very similar heliocen- 
tric velocity). The two estimates are in excellent agreement, which 
is a robust check on the reality of the signal. Fig. [2] also shows 
la error bars for the estimated power spectrum. The errors were 
computed by including both the effect of having a fine number of 
estimates of the true spectrum ("cosmic variance") as well as hav- 
ing noisy visibility measurements. The first term was estimated as 
Phi (U) \/D/2itU, while the second term was taken from the mea- 
sured fluctuations of Pm(U) for channels with no HI emission. 
The total la error bars were determined by combining the two con- 
tributions in quadrature. At small baselines (U ^ 5kA), the first 
term dominates the error budget, while the noise on the measured 
visibilities dominates at large baselines. 

If the observed power spectrum is due to turbulence in the in- 
terstellar medium, the slope of the derived power spectr um could 
chang e depending on the width of the frequency channel iLazarianl 
( 1995)). Turbulence gives rise to fluctuations in both the density 
and velocity of the gas and both of these contribute to the observed 
intensity fluctuations. Lazarian & Pogosyan (2000) show that the 
statistical properties of the velocity and the density fluctuations 
can be disentangled if the slope of the observed power spectrum 
changes depending on whether the intensity is averaged over a ve- 
locity range that is large ("thick slices") or small ("thin slices") 
compared to the turbulent velocity dispersion. For DDO 210 the 
observed velocity dispersion is ~ 6.5 km s _1 and is approximately 
constant across the galaxy. Pm(U) estimated from visibilities av- 
eraged across the 10 channels over which the emission has a simi- 
lar spatial distribution is shown in Fig.|5] Within our measurement 
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uvdistance (Kk ) 



Figure 2. Absolute value of the real part of Pui(U) plotted as a function 
of uv-distance U for the VLA data (dash line) and the GMRT data (solid 
line) for a frequency channel with significant HI emission. The GMRT data 
is exactly the same as in Fig. also the VLA frequency channel has the 
same heliocentric velocities as the GMRT data. 
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Figure 3. Absolute value of the real part of Pjji (U) determined after aver- 
aging the GMRT visibilities over the 10 central frequency channels which 
contain significant HI emission. The best fit power law using the UV range 
2 — 12 kA is shown as a solid line. The best fit convolved power law over 
the UV range — 12 kA is shown as a dotted line. See the text for more 
details. 



errors, we find no change in the power spectrum with increasing 
velocity width in the range from ~ 1.7 km s _1 to ~ 17 km s . 
We discuss the implications of this in Sec|5| but we note here that 
in this means that one can work with the emission averaged over 
channels to improve the signal to noise ratio. 

The power spectrum Pm(U) shown in Fig. [5] appears to fall 
like a power law beyond a uv distance of ~ 1 kA. The best fit power 
law was determined by minimising 



X 



E 



pm(Ui) - AU?\ 



(11) 



with respect to the amplitude A and the power law index a. It 
is also necessary to fix the uv range while doing the fit. As can 
be seen in Pig. [3] Pm(U) flattens below 1 kA; at these scales 
one would expect that the convolution with the window function 
|Wi/(t/)| 2 would be important. If we assume that the convolution 
with the window function can be ignored for U > 1 kA and do 
the minimisation described in Eqn. 1111 with the lower limit fixed 
at £/min = 1 kA, but several different upper limits (7 ma x, we find 
that the minimum x 2 /DOF decreases up to t/ max = 12 kA but not 
beyond. For the range 1 — 12 kA the best fit power law index is 
q=— 3.55lo'3Q. Similarly, keeping U m in = 2 kA, the best fit power 
law index for the uv-range 2 — 12 kA is a=— 2.85^o'g5- The best 
fit for the uv-range 2 — 12 kA is shown in Fig. [3] We note that the 
best fit a is largely determined by the smallest baselines U fa Umin 
as these have the smallest error bars, and it does not change signif- 
icantly if (7 m ax is changed to say 6 kA or even 50 kA. As one goes 
to larger and larger baselines however, the minimum in the % 2 dis- 
tribution gets more and more shallow - equivalently the error bars 
on a get larger and larger. 

To check the assumption that the convolution with the win- 
dow function can be ignored for U > 1 kA, we modelled the 
window function using a Gaussian W v 0) = e~ u /2D of width 
D = 0.8 kA. For a power law form of Phi (U) we find that the con- 
volution actually produces measurable deviations from the power 
law at U ^ 2 kA. The observed power spectrum is steeper than 
the underlying power law spectrum in the range 1 — 2 kA, while at 
baselines U ^ 1 kA it is flatter than the underlying power law. We 
hence repeated the least squares fit to the observed power spectrum 



using a power law convolved with the assumed Gaussian window 
function as the input model. For this model the best fit slope for 
the range U < 12kA is a = -2.1+ ;° . The value of the best 
fit slope when including baselines < 2kA could depend on how 
well the true window function has been modelled. In our analysis 
we have approximated this with a Gaussian; any deviations from 
a Gaussian could introduce unknown biases in our estimate of the 
slope. With this consideration in mind, we also did fits for the re- 
stricted uv-range 2 — 12kA. For this uv-range, the best fit slope is 
a — —2.75 ± 0.45. While the fit is good over the fitting range, the 
model substantially over estimates the power on large scales. This 
is consistent with the flatter slope (~ —2.1) obtained when fitting 
over the entire U ^ 12kA range, and might be indicative of a devi- 
ation from a power law on large scales. However, given the uncer- 
tainty in the modelling of the window function it is difficult to make 
a quantitative statement. On the other hand, our modelling does in- 
dicate that beyond about ~ 2kA the effect of the window function is 
marginal. The value of the slope obtained after taking into account 
the convolution (viz. —2.75 ± 0.45) is consistent within the error 
bars with the value a — — 2.85^0 65 obtained without accounting 
for the convolution. We hence adopt the value a = —2.75 ± 0.45 
as slope of the intrinsic power law over these spatial scales. 



5 DISCUSSION 

Over the uv-range 2 — 12 kA (corresponding to spatial scales of 
80pc - 500pc) the power spectrum of HI intensity fluctuations in 
the dwarf galaxy DDO 210 is well described by a power law with 
slope a = —2.75 ± 0.45. This value is very close to the slope 
of ~ — 3 estimated from the HI emission in our own galaxy (on 
spatial scales of ~ 5 pc to 200 pc) by Crovisier & Dickey(1983) 
and Green(1993), and for the SMC (on spatial scales between 
~30 pc to 4 kpc; Stanimirovic et al. 1999). The reason why the 
HI emission power spectrum has a similar slope in the Milky- 
way as in a galaxy with more than 3 orders of magnitude less 
HI mass suggests that some fundamental physical process is re- 
sponsible for determining the power spectrum slope. Power spec- 
tral slopes ~ — 3 are expected from optically thick turbulent gas 
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lLazarian & Pogosvanl l2004h ). The HI emission from late type 
dwarfs like DD O 210 is however expected to be opticall y thin, 
(see e.g. lHavnes & Giovanni 1984hlGiovanelli et alll9 941. Mod- 
els like those of Goldman 1 200G) which give a spectral slope ~ — 3 
in optically thin gas may be more relevant in this case. 

Unlike the Milky way, DDO 21 has no spiral arms and it also 
has no measurable star formation I van Zee (2000)). The lack of 
dependence of the power spectrum slope on the star formation rate 
that we find is similar to the finding of Dib et al. (2006) that the 
characteristic HI velocity dispersion in galaxies shows a minimum 
level that is independent of the star formation rate. In this aspect, 
the turbulence in the gas distribution appears to be different from 
that in the stellar distribution, where the power spectrum's slope 
does seem to depend on the star formation rate (Willett 2005). It is 
possible that the energy source driving the HI turbulence is gravi- 
tational energy input on large scales. Detailed simulations suggest 
that large scale gravitational energy input could explain the sub- 
stantial turbulent velocity dispersion seen in the outer parts of the 
extended gas disk of NGC 2915 (Wada et al. (2002)) as well as 
for the rich structures seen in Holmberg II (Dib & Burket (2005)). 
Gravitational energy input on large scales (perhaps by tidal torques 
from other local group members) would also be consistent with the 
asymmetric large scale HI distribution in DDO 210. 

As mentioned in Sec.|4| turbulence produces fluctuations in 
both the density and velocity fields, in turn, both of these contribute 
to fluctuations of the HI intensity. Lazarian &Pogosyan (2000) 
show that the slope of the HI intensity power spectrum goes like 
n + C2/2 when measured for channel widths that are "thin" com- 
pared to the turbulent velocity dispersion, while for channel widths 
that are "thick" compared to the velocity dispersion the HI intensity 
power spectrum's slope is n, where n is the slope of the power spec- 
trum of density fluctuations and £2 is the spectral slope of the ve- 
locity structure function. Essentially, for velocity widths large com- 
pared to the velocity dispersion ("thick" slices) the velocity infor- 
mation gets averaged out and the intensity fluctuations represents 
density fluctuations. Recall that we find that within the measure- 
ment errors the slope of the HI power spectrum does not change 
for velocity widths ranging from ~ 1.7 km s _1 to ~ 17 km s _1 . 
The observed HI velocity dispersion of 6.5 km s _1 has contribu- 
tions from both thermal and turbulent random motions. Our "thick" 
channel is wider than the total velocity dispersion, and hence we 
can be sure that the detected power spectrum represents the power 
spectrum of the HI density fluctuations. Since we get the same 
slope ~ —2.75 ± 0.45 even for ~ 1.7 km s _1 channels, either 
the turbulence velocity dispersion is smaller than 1.7 km s _1 or 
£2 ~ ± 0.9. Note that the latter possibility is consistent with 
£2 = 2/3 which is predicted by the Kolmogorov theory of turbu- 
lence. 
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